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Abstract 



We study the migration of chemotactic wild type Escherichia coli popu- 
lations in semi-solid ('soft') agar in the concentration range C = 0.15 — 
0.5%(w/v). For C < 0.35%, expanding bacterial colonies display character- 
istic chemotactic 'rings'. At C = 0.35%, however, bacteria migrate as broad 
circular bands rather than sharp rings. These are growth/diffusion waves 
arising because of suppression of chemotaxis by the agar and have not been 
previously reported experimentally. For C = 0.4 — 0.5% expanding colonies 
do not span the depth of the agar and develop pronounced front instabili- 
ties. The migration front speed is weakly dependent of agar concentration 
below C = 0.25%, but decreases sharply above this value. We discuss these 
observations in terms of an extended Keller-Segel model for which we de- 
rived novel transport parameter expressions accounting for perturbations of 
the chemotactic response by collisions with the agar. The model allows to 
fit the observed front speed decay in the range C = 0.15 — 0.35%, and its 
solutions qualitatively reproduce the observed transition from chemotactic 
to growth/diffusion bands. We discuss the implications of our results for the 
study of bacteria in porous media and for the design of improved bacterio- 
logical chemotaxis assays. 

Key words: Escherichia coli; motility; semi-solid agar; porous media; 
chemotaxis; population dynamics 
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Introduction 

Much is understood about the motihty of flagehated bacteria in open hq- 
uid media (jl|) and on sohd surfaces ([3). In contrast, bacterial locomotion 
within semi-solid media is much less well studied, even though bacteria of- 
ten colonise three dimensional semi-solid environments, e.g. host tissues or 
foodstuffs. 'Soft' agar with concentration 0.2% < C < 0.35% (through- 
out, % = %w/v), a gel network whose main component is the semi- flexible 
polysaccharide agarose ([3), is a plausible model for many of these kinds 
of solid environments. It has been used in studies aimed at understanding 
motile microbial pathogens growing inside the semi-solid matrix of a vari- 
ety of foods and of infected hosts, and therefore for predicting spoilage and 
infection (0-0). 

Soft agar was first introduced into microbiology for assaying chemotaxis 
(0). Chemotactic wild-type Escherichia coli inoculated at one end of a 
capillary filled with nutrient buffer spread out in 'bands' as they succes- 
sively deplete the medium of oxygen and various nutrients. In a Petri dish 
filled with nutrient agar, the observation of successive sharp circular bands 
('rings') progressing outwards from the colony inoculated into the centre 
of the soft agar is taken to confirm the chemotaxis genotype, since non- 
chemotactic mutants spread out uniformly (j^). 

Interestingly, the agar concentration in this widely-used chemotaxis as- 
say is not standardized, varying from investigator to investigator, or even 
within the same study, in the range of 0.2% < C < 0.4%. The tacit as- 
sumption seems to be that, as long as concentrations are in the 'soft range', 
agar conveniently suppresses thermal and biological convection in the liquid 
medium which hosts bacteria, but otherwise has no interesting effect. Wolfe 
and Berg's investigation (0) of a number chemotactic mutants in soft agar 
(C = 0.2% — 0.35%) appears to confirm this assumption. They report no 
qualitative diff'erence except a slowing down of the spreading front. 

However, the run-and-tumble motion of E. coli ((ll) and similar bacteria 
suggests a priori that soft agar should affect chemotaxis. The pore size 
of soft agar is ~ 1/um (fiol ). the same order as typical bacterial run lengths. 
Since cells perform chemotaxis by altering their tumble frequency and there- 
fore run length, the structure of agar may therefore interfere with the ability 
to chemotactically 'bias' random walks. The expectation is that this inter- 
ference will be more pronounced in agar than in other porous media, such 
as sand and soil, where most pore sizes (~ lO^m to 1 mm) are very much 
larger than typical run lengths. Indeed, a study of Pseudomonas putida in 
sand columns with grain sizes in the range 80 — 800/xm found no effect on 
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chemotaxis (j^. 

In this work, we show by experiment and theory that the chemotaxis 
of E. coli in soft agar is indeed strongly affected as the concentration is in- 
creased in the range C = 0.15 — 0.5%. We observe that, as C is increased, 
the colony growth pattern changes qualitatively. Also, the speed of the mi- 
grating front is weakly dependent on C at C < 0.25%, but decreases sharply 
with agar concentration above this value. To understand our observations. 



we formulate a modified Keller-Segel type model (jlll . Il2l ) with transport 



parameters that are functions of agar concentration. These functions are 
derived by extending a recent description of bacterial chemotactic response 

0), 

and account for the diminished ability of bacteria to detect chemical 
gradients in dense gels. Our model is able to fit quantitatively the observed 
dependence of front speed on agar concentration. 



Methods 

We used the E. coli K-12 derivative AB1157, a chemotactic wild-type 

0). 

Plate cultures from frozen stocks were refrigerated at 4°C for up to 3 weeks 
prior to use. Luria Broth (LB) agar was prepared by adding 1.5-5 g/1 of 
Difco Bacto-agar to LB BrothflO g/1 Difco Bacto-Tryptone, 5 g/1 Difco 
Yeast extract, 10 g/1 NaCl) (jlSl). The mixture was autoclaved at 121°C 
under 1.02 atm for 30 min, and left to cool for an hour. Agar plates were 
then prepared by pouring hot (45-50° C) nutrient agar into standard sized 
(100 mm diameter, 12 mm deep) plastic Petri dishes (Sterilin) on a level 
surface; 58 ml were poured in each plate (final agar thickness: lO.OitO.l mm). 
Poured plates were left to set for about a day at ambient temperature (22- 
27°C). The final pH of the agar was 7.3 ± 0.1. 

Late exponential phase cultures (ODgoo = 0.8 — 1.3) were prepared by 
inoculating single colonies in LB-filled flasks incubated at 30.0°C and shaken 
at 200 rpm. Cultures were then diluted to ODgoo = 0.1 (r^ 10® viable 
cells/ml) and agar plates were inoculated by delivering a droplet of < Ifil 
via a 2/il pipette. Inocula were left to sit on the agar for about an hour; 
then a thin layer (thickness ~ 1.5 mm, corresponding to ~ 9 ml) of sterile 
filtered mineral oil (Sigma) was poured on the plates. Except at the highest 
agar concentrations studied (C > 0.4%), the inoculum had been (to the 
naked eye) completely assimilated by the time the oil was poured, with no 
visible evidence of the pouring spreading organisms to other parts of the agar 
surface. At the highest C, the inoculum was not significantly incorporated 
into the agar after an hour, so that careful pouring was needed to minimise 
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spreading. When significant spreading did occur, the plates were discarded. 
The thin oil layer kept evaporative losses to < 1% (weight) during our 
observations, but did not generate anaerobic conditions (|l6l ). 

The Petri dishes were incubated at 30.0±0.5°C on a dark background, il- 
luminated from the sides, and imaged at 30 minute intervals using computer- 
controlled CCD cameras. Images were analysed using ImageJ (NIH) and 
IDL (RSI, Boulder, CO): a, static noise (by subtracting the first frame in 
each sequence) was removed; b, slightly non-uniform illumination was cor- 
rected for (subtracting the image background fitted using a sliding paraboloid 
with 20 pixels rolling ball radius); c, histograms were matched histograms 
using gray-scale mapping (llTI ). Then we obtained azimuthally-averaged ra- 
dial intensity profiles from those images without significant blebs (see below 
for observations and further discussion on blebs). The images presented in 
Fig. [1^ and Fig. [2] were processed using only steps b and c. 

At high C, image thresholding enabled us to determine the colony area, 
A, from which we calculated the radius r = A/tt. At lower C we fitted a 
circle to the intensity maximum in each image and determined the area of 
the fitted circle, from which r was then calculated. 



Model of chemotactic E. coli populations in agar 

The fundamental processes in agar plates inoculated with bacteria are growth 
due to nutrient uptake and dispersion due to chemotactic motility, which can 
be modelled by generalized Keller-Segel models (I 111 . Migrating popu- 



lations of bacteria in agar have been described for bacteria chemotactically 



sensing nutrients ([l9|, l2d ) or attractants secreted by the cells themselves 



(|2ll). However, these models ignore, or inconsistently account for, the effect 



of agar on bacterial chemotaxis. Since agar is a porous gel, one might think 



that existing descriptions of other porous media ([23, |23| ) should be appli- 
cable to agar. We will show below that these descriptions are incorrect. 
The model presented here is an adaptation of that originally formulated by 
Lauffenburger, Kennedy and Aris (LKA) for a chemotactic population with 
growth in one dimension (fl^ ). Our model differs from the LKA model in 
three ways. First, we model growth as logistic, while LKA used a linear 
term. Secondly, we work in two dimensions, since LKA's ID analysis is not 
adequate for modelling Petri dishes for early times. Finally and most cru- 
cially, bacterial transport coefficients in our model are not constants, but 
functions of agar concentration derived from a recent model of chemotactic 
response (fl^). 
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Model equations 



The starting point equations of our model are: 



db 
di 
ds 
di 



-V- 



[-^(s, C)V6 + v^(., C)b] + k,b (^g{s) - ^) (1) 
5 - ^kgg{s)b. (2) 



Equation [T] expresses the conservation of bacteria, with population density 
b{r,t). This population evolves in response to the combined effect of its 
diffusive and chemotactic fluxes, with diffusivity fi{s, C) and drift velocity 
v-^(s,C); these are in general functions of both substrate and agar concen- 
trations, s(r, t) and C respectively. The bacterial population also evolves by 
growth, with birth rate kgg{s), where kg is the maximum growth rate and 
g{s) is a function of substrate uptake, and a death rate —kgb/kf,, where ki, 
is the carrying capacity of the population. Equation ([2|) models the conser- 
vation of the first, most readily metabolised substrate, with concentration 
s{r,t), diffusing with diffusivity D and being consumed by bacteria at a 
rate —kgg{s). In tryptone broth or LB this substrate is L-serine (0, I20I ). 
In the consumption term of ([2]), y is the bacterial yield upon consump- 
tion {b = Ys). We now assume the following relations: iJ,{s,C) = fJ,{C); 



v^(s,C) = x{C)yfx{s), where /^(s) = ^q^; and g{s) = That is, we 



assume the diffusivity to be isotropic and independent of substrate concen- 
tration; n depends only on the spatially uniform agar concentration C (see 
below). The chemotactic velocity, v^, is assumed linear in the gradient of a 
'receptor-adsorption' function fxi^)', ky. is the characteristic saturation con- 
centration of the chemotactic response (0). The proportionality constant 
is the chemotactic coefficient x which, like ^, is assumed to depend only on 
agar concentration. The relations involving and are approximate forms 
valid in the limit of shallow concentration gradients (jisi ). Bacterial growth 
depends on substrate concentration through a Monod-type growth function 
g{s) (0); kg is the characteristic saturation concentration for growth. 

We have derived the dependence of the diffusivity, //, and chemotactic 
parameter, on agar concentration by modifying de Gennes' integral model 
of bacterial chemotactic response (|l3l ). as detailed in Appendix [XI The 
model quantifies the intuition that in a dense matrix of obstacles bacteria 
are reoriented by collisions with the matrix as well as by tumbles, making 
chemotaxis inefficient. Such collisions increase the effective bacterial tumble 
rate from the in liquido value oq to a{C) = ao(l + f{C), where the function 
/(C) quantifies the collision-induced concentration-dependent increase of 
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the tumble rate (see below). From our model it can be shown (see Appendix 
E| that the chemotactic transport parameters in agar are given by: 

/i(C) = /.o [1 + f{C)]-' ; x{C) = XO [1 + f{C)]-' [/(C)] (3) 

where and xo are the bacterial diffusivity and chemotactic coefficient in 
the absence of agar. The agar concentration-dependent integral is given 
by equation ()A16p and its value depends on the form of the chemotactic 
response function, K(t), see equation ()A17p . The function /(C) gauges the 
increase with agar concentration of the tumble rate in agar, a^(C), with 
respect to its in liquido value uq. That is, we assume a^(C) = aof{C). Since 
collisions are more frequent for a higher density of obstacles, we expect /(C) 
to monotonically increase with C. We adopt the ansatz /(C) = exp[ '^^^^ ], 
where Co is a characteristic concentration. The concentration 'shift' Ci 
accounts for the possibility that the tumble rate in agar can recover its in 
liquido value for small, but nonzero agar concentrations: a^iC < Ci) — )■ 
so a{C < Ci) ao, Ci > 0. 

To further understand the experimentally observed migration transition 
and to compare our results to those derived for bacteria in porous media 



22 . |23| ). we have also derived asymptotic limits to the expression ([3]) for 



x(C). As shown in Appendix El these asymptotic limits are: 

/xo[l + /(C)]-2[l-K/(C)] ifd(C)«l 
1/3X0 [l + /(C)]-=^ ifd(C)»l. 

where k and /3 are constants (see Appendix E]) and a(C) = 1 + a^/ao = 
1 + /(C) is the dimensionless effective tumble rate in agar (see (|A6P ). The 
limits dH) reflect the effect of confinement in agar on chemotaxis. At low 
concentrations, agar does not significantly impede chemotaxis, and bacte- 
ria can tumble relatively freely: a{C) ~ 1 ('efficient' limit). At higher 
concentrations, frequent collisions with the agar d(C) ^ 1 confuse the 
chemotactic response (the 'confused limit'). Neither expressions ([3]), nor 
the limits coincide with those derived in previous models of bacteria in 
porous media. These models treat bacteria in porous media like molecu- 
lar gases, and so derive transpo rt p arameters obeying the balance relation: 



MCO _ xM - [1 + dad, 1231). However, in our model this rela- 



tion does not hold (even for very low C) because of the effect of collisions 
on chemotaxis. It is only recovered in the absence of agar. We assume 
changes in the swimming speed, D or kg with agar concentration are negli- 



gible (1271, 12 
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Geometry, scaling and model parameters 

We will consider only the two-dimensional, axisymmetric limit of our equa- 
tions. Our assumption is that, modulo a time shift, fully developed bacte- 
rial front dynamics are insensitive to particular initial conditions (provided 
the initial colony is azimuthally symmetric). The characteristic length and 
time scales for our experiments are millimeters and hours, so we rescale our 
equations by Tg ~ k~^, where kg is the growth rate, and Ig y^^o/kg, 
the average length a cell diffuses during a doubling time (in the absence of 
agar). The population density is rescaled by its carrying capacity, k^,, and 
all concentrations by the initial substrate concentration, sq. We rescale all 
diffusivities by that of the bacterial population (in the absence of agar), fiQ. 
We also rescale the yield Y by the maximum possible yield k^/sQ. Thus: 
r = r^kg/iiQ] T = kgt; B = h/h; S = s/sq; N = D/hq] Sq = Xo//^o; = 
kjso; = ks/so; H = h/{Y so); M{C) = //(C)///o; X{C) = x{C)/xo- 
The model equations ([TJ and ([2]) in dimensionless form then read 

^ = M{C)V^B - 5oX{C)V ■ ^V5) + B [G{S) - B] , (5) 
^ = NV^S - HG{S)B, (6) 

where F^{S) = sTk^,G{S) 

— s+Ks ^^'^ where the dependence on agar 
concentration is through the functions: 

M(C) = [1 + fiC)]-' ; XiC) = [1 + fiC)]-' [/(C)] , (7) 

where we recall /(C) = exp[^^^^p-], where Co and Ci are the characteristic 
concentrations introduced earlier. The parameter 5q is significant: it mea- 
sures the relative magnitude of chemotactic advection to random diffusion 
in the absence of agar (a chemotactic Peclet number). It is how this ratio 
is modified by agar which leads to surprising results, as we will see. We 
have also defined the dimensionless parameters N, the ratio of nutrient and 
bacterial diffusivities, and H, the ratio of carrying capacitance of the bac- 
terial population to the maximum population obtainable from the nutrient 
available. Equations ([5]) and ([6]) are subject to no-flux boundary conditions 
and to the initial conditions: 

B{R,0) = e~^] S{R,0) = l-e~^, (8) 

where a is the width of an initial Gaussian packet of bacteria. The fol- 
lowing parameter values were used to solve our equations: growth rate, 
kg = 0.7h~^ (from in liquido growth curve); initial cell concentration, 6o(= 
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/cfe) = 3.5 X 10^ cells/ml (from viable counts); initial substrate concentration 
(of L-serine in LB), sq = ImM (5-8 mM (|29l )): cell diffusivity (no agar), 
/io = 5.7mm^/h (1.2mm Vh B); chemotactic parameter (no agar), xo = 
600mm^/li (450mm^/h to a-methylaspartate (|3d ) ) : substrate diffusivity, 
D = 3mm^/h; chemotactic threshold concentration, = 0.5 mM (0.2 mM 
for a-methylaspartate (|3ol)): growth threshold concentration, ks = ImM; 
yield, Y = 10^^ cells/ml/M. The growth rate and initial cell concentration 
were determined by our own independent experiments indicated in brackets. 
All other parameters are based on experimental literature values for E. coli, 
many of which have been used in other models of E. coli migration (|3ll . |32|). 
Reference literature values close to parameters we changed are reported in 
brackets above. In addition to these macroscopic parameters, we use the 
in liquido tumble rate ao = 1 s and the constant Aq = 0.5 to calculate the 
integral in ^ using ()A16P and ()A17p (jsil ). The concentrations Co and 
Ci are free parameters, fixed by fitting the predicted front speeds with those 
we observed experimentally (see Results and Discussion). Prior to perform- 
ing the fit, the values of the parameters kg, fiQ, xo '^sre adjusted 
slightly to match the values of experimental and predicted band speed for 
C = 0.15% (assuming this is the same as in liquido). 

With these parameters, the dimensionless constants of the model have 
the values: 6o = 105; = 0.53; Ks = I; N = 0.5; H = 3.5. The above 
parameters will not be changed in our investigation and the initial packet 
width a is fixed at 2. Equations (jS]) and ^ in ID axisymmetric form were 
solved numerically for C = 0.15 — 0.35% on a linear domain (L = 100) using 
Matlab subject to initial conditions dH) and no-flux boundary conditions. 
Migration front speeds were obtained by subtracting the position of the 
leading edge inflexion points of solution proflles calculated at neighbouring 
time points and dividing by the time interval. Like in experiment, these 
speeds were calculated in the linear growth regime (long times) where speed 
does not change with radius. 



Results 

Observations on migration morphology and radial dynamics 

We first report qualitative features of colony morphology and dynamics. 
For all concentrations studied (C = 0.15 — 0.5%) it takes 5-7 hours for the 
bacterial inoculum on the agar surface to become visible. The inoculum 
then grows in optical density and, after an additional time lag of 1-50 hours 
(likely caused by the oil overlay, but with no infiuence on the reproducibility 
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of subsequent front dynamics), the initial bacterial colony migrates across 
the plates. Stills from early and advanced stages of colony migration for 
concentrations in the range 0.15% < C < 0.35%, are shown in Fig. [1^. 
Two striking effects of increasing concentration are immediately apparent: 
the change from a morphology displaying sharp rings to one which is more 
diffuse and featureless, and the loss of circular symmetry in the advanced 
stage of migration at high concentrations {C f« 0.35%). 

At the lowest concentrations sampled, C = 0.15 — 0.2%, bacteria migrate 
as sharp circular bands inside the agar. We observed two bands in succession. 
The first band sharpens as it migrates across the plate. Fig. [TJ),c. The 
second band is slower than the first and also appears to sharpen as it travels. 
Interestingly, the first band at C = 0.15% initially displays internal structure 
(a double band, see first frame of Fig. [lb) and is reflected from the plate 
walls (not shown) before the second band catches up with it. 

Bacteria also migrate as circular bands for C = 0.25 — 0.3%. Again 
two bands were observed, but now they travelled together Fig. [T]l,e. At 
C = 0.35%, sharp bands are no longer visible (Fig. [T]F). The colony grows 
from the inoculum as a circular disk with a slightly nebulous front (Fig. [1| . 
The intensity across the disk is initially approximately uniform, falling off 
at the edges, defining the band front Fig. [T]F. At later times, however, it 
displays a broad band structure. We did not follow the radial development 
of these bands to the edge of the plate because the colony front develops 
instabilities (blebs) that disrupt circular symmetry. 

Visual inspection confirmed that bacteria had spread from the surface 
inoculum into the agar to a significant depth. For concentrations supporting 
bands these are initially hard to resolve for radii smaller than the agar depth: 
the colony appears like a uniform expanding circle from above (Fig. [1^, top 
row). For larger radii the first band is visible and clearly spans the depth of 
the agar, as observed by Adler (0). For colonies with two distinct bands, it 
is not clear at what depths the second bands occur; from our images they 
seem to be further inside the agar. Microscopy (not shown) reveals that for 
C < 0.4% bacteria penetrate significantly beyond 1mm in depth, but for 
this and larger concentrations it seems that agar limits penetration to a few 
mm from the surface. 

At C = 0.4 — 0.5%, shown in Fig. [21 the expanding colony appeared 
as homogeneous, solid circles initially. However, the front invariably de- 
veloped extensive blebs. The blebbing instability set in earlier for higher 
concentrations (e.g. at C = 0.4% blebs appeared when the colony radius 
was beyond a third of the plate radius, while at C = 0.45% it appeared at 
around one sixth). The blebs developed into wedge shaped sectors, giving 



E. coli in soft agar 



10 



the colony an overall flower shape (Fig. [2]). At these concentrations, the 
colony also appeared to spread on the surface of the agar (though not by 
classical 'swarming'), but we did not investigate such surface migration. 



Effect of concentration on radial migration 

In Fig. [3] we plot the radius of the outermost migrating front (band), r, 
against the time, t = ti — At;, elapsed since inception of visible colony 
growth, where ti is the time since inoculation and At; is the latency time 
before a colony grows out. We estimated At; from the intersection of a linear 
fit to the raw radial data with the time axis (Fig. [3l inset). A substantial 
portion of the radial growth is linear in time for 0.15% < C < 0.4%. Linear 
portions can also be identified for 0.45 and 0.5% (not shown), though the 
extent of these data is severely limited by the formation of blebs. Slopes from 
the fits to the radial growth curves in the range 0.15 — 0.35% (0.15 — 0.5%) 
are plotted as a function of agar concentration in Fig. |4] (and inset). At 
C < 0.25% the migration speed is at best weakly affected by concentration, 
but beyond this value it decreases dramatically. Our model can account for 
this behaviour (see below). 



Theoretical front speed decay and band profiles 

A fit to the experimental front speed data from solutions of our full model 
using relations ^ is shown in Fig. HJ also shown are the 'efficient' and 'con- 
fused' limits of the model for the same parameters. As expected, the efficient 
(confused) limit is a good description at low (high) concentration. The evo- 
lution of the theoretical band profiles corresponding to the full model best fit 
is shown in Fig. [5] (left). Also shown in Fig. [5] (right) is the prediction using 
transport parameters from gas kinetic models derived for bacteria in porous 
media: = = [1 + (22, 23). As concentration is increased 



in the experimental range C = 0.15 — 0.35%, the full model band profiles 
displays a gradual transition from sharp, chemotaxis-dominated bands to 
broader, growth/diffusion-dominated bands. In the gas kinetic model, be- 
cause the chemotaxis parameter, ^^id diffusivity, /i, have the same func- 
tional dependence on C, profiles remain sharp for all concentrations. The 
rounded profiles predicted by our model arise from suppression ('confusion') 
of the chemotactic response caused by bacterial collisions with the agar. 
When the chemotactic flux in equation ([1]) becomes negligible with respect 
to the fluxes due to logistic growth and diffusion, the travelling band solu- 
tions to ([1]) and ([2]) change from sharp, fast chemotaxis-dominated bands to 
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slower, broader bands driven by growth/diffusion processes. This is what we 
observe experimentahy. The breakdown of the model in the range C = 0.4- 
0.5%, evident from the inset to Fig. U is explained in the discussion below. 



Discussion 

We have experimentally studied the migration of chemotactic E. coli popula- 
tions in soft agar of concentration in the range C = 0.15—0.5%. Consistently 
with other investigators we find that increasing agar concentration decreases 



the speed of propagation of the bacterial front ([a, |a, |3J) and severely ham- 
pers penetration for C > 0. 5% (0, S). However, our work also reveals a 
hitherto unobserved transition in the dynamics of the population as agar 
concentration increases. The gradual transition is from a dynamics display- 
ing characteristic sharp chemotactic bands (rings) to one where the bacteria 
travel as broader bands. By increasing the chemotaxis to diffusion ratio 
^0 = Xo/fJ'O Lauffenburger et al. (LKA) theoretically studied the transi- 
tion from sharp chemotactic to broader growth/diffusion bands, but failed 
to find evidence for the latter in studies of chemotaxis in capillaries 0). 
Interestingly, we have discovered that sufficiently concentrated agar pro- 
vides an environment where chemotaxis is suppressed and growth/diffusion 
processes can be observed to dominate the band dynamics. 

To understand our experimental results we also built a model of bacte- 
rial migration in agar. We extended the LKA model and coupled it to the 
first full expressions for the concentration dependence of bacterial diffusiv- 
ity n{C) and chemotactic parameter x{C) in agar. We derived these (see 
Appendix |A]) from an adaptation to agar of de Gennes' model of bacterial 
chemotactic response 

0). 

Collisions with the matrix of concentrated agar 
(effective tumble rate a{C) = ao(l + f{C)), where /(C) = ex.p[^^^f^]) con- 
fuse this response causing fJ,{C) and xiC) to have different functional forms. 
Our model can thus predict the band transition we observe experimentally. 
We obtained a best fit of the theoretical front speeds to the experimentally 
observed values (Fig. H]) in the concentration range 0.15-0.35%, finding the 
characteristic concentrations Co = 0.035% and Ci = 0.28%. In comparing 
model profiles. Fig. [5l with experimental ones, Fig. [TJ)-f, we note that the 
vertical axes in the latter probably do not map linearly to cell density due to 
multiple scattering effects. In addition, 'dead' or non-motile bacteria con- 
tribute to the signal but do not contribute in theoretical plots. With these 
caveats, we see that for C = 0.15-0.35% our model qualitatively reproduces 
the experimentally observed transition in the colony (band) profile at long 
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times rather nicely (Fig. [5l left): bands change from sharp to broad as 
concentration in increased. 

The model breaks down for C = 0.4-0.5%. At these concentrations bac- 
terial diffusivity becomes very small (e.g. M(0.4%) = 0.03) and equations 
([TJ and ([2]) predict a front speed independent of C. However, the measured 
(early, bleb-less) front speed continues to fall sharply with C, see inset to 
Fig. m One reason the model fails is that small diffusivity affects growth at 
high C. During a doubling time bacteria in 0.4% agar are able to diffuse ~ 6 
times less far than for C = 0.15% {in liquido), which increases competition 
for nutrients with neighbours. Further, at high C small bacterial diffusivity 
means growth is limited by that of nutrients: /Lt(0.4%)/L' = 0.06. Diffusion 
limited grow th is known to produce branching instabilities like those we 
observe (|35l ).To fully explain high concentration colony morphologies (Fig. 
[2|) changes in gene expression in response to low nutrient levels will also 
need to be considered. An interesting possibility is that in high C agar cell 
densities could reach large enough values to elicit quorum sensing responses 
(js^). Experimentally, the situation for C > 0.4% is also complicated by 
the observation of coexisting subpopulations (see results and also (0)), one 
growing on the surface and one in the bulk, which does not penetrate very 
deeply (the dynamics is no longer 2D as assumed). Modelling these very 
different conditions is left to a future study. 

We have so far been implicitly discussing the first (front) band. Ex- 
perimentally a second band is also observed for C < 0.35% which, as agar 
concentration is increased, travels closer and closer to the first (see Fig. 
[TJj-f). As mentioned bacteria in LB preferentially metabolise one nutri- 
ent at a time: the first band aerobically consumes L-serine and the second 
L-aspartate, with a roughly constant metabolic delay Tm between bands 
(0, S, [20I ). Thus the maximum spacing between bands Li, ~ VF{C)Tm will 
decay with C like vf{C), the speed of the first band. In this paper the em- 
phasis has been on explaining the experimentally observed shape transition 
of the first band. In the future, it will be interesting to extend our model 
and experiments to quantify chemotaxis and its suppression for all nutri- 
ents consumed. Accounting for multiple bands, as well as using improved 
receptor-adsorption functions for growth and chemotaxis, will allow more 
realistic predictions for the trailing edge of the bands. 

The suppression of chemotaxis we have studied is relevant to the migra- 
tion of bacteria in porous materials other than agar, important in bioreme- 
diation (22) and food spoilage (Q). As discussed above, previous gas kinetic 
models of bacterial migration in porous media (2j, |23| ) do not account for 
the possibility of the chemotactic response becoming confused by collisions 
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with agar. This neglect, which is an imphcit consequence of assuming bac- 
terial populations behave like molecular gases, invalidates the predictions of 
these models in porous media with a finite concentration of obstacles, even 
if dilute. Gas kinetic models can provide good fits to our experimental front 
speed data (with different values for the characteristic concentrations Cq 
and Ci), but cannot also reproduce the experimentally observed transition 
in front shapes. On the other hand, provided pores are larger than a cell, our 
model accurately describes the transport of chemotactic bacteria in general 
porous media. 

Our results also have potentially important implications for microbio- 
logical practice. Microbiologists studying motility often make chemotactic 
mutants which are screened for using chemotaxis assays. One of these assays, 
the 'motility assay', involves inoculating soft agar and imaging the resulting 
bacterial colony, like we have done in this study. The agar concentration 
for such assays is not standard (values in the range 0.1-0.4% can be found 
in the literature (j37l . ISg)). and seems to be a matter of convenience (e.g. 
larg er concentration allows to study more than one colony in the same plate 
(j39l )). When chemotactic mutants are screened for, the chemotactic band 
phenotype is sought for as a marker of chemotaxis, its absence denoting a 
successful chemotactic mutant (0) or a failed restoration of the chemotaxis 
phenotype (js?! ). Our experiments suggest, however, that chemotactic run- 
and-tumble bacteria above a certain (still soft) concentration of agar will 
fail to show the band phenotype. Thus, if agar plates are used to assay 
for chemotaxis it will be important take into account the possibility that 
suppression of the band phenotype by the physical environment may oc- 
cur. Performing assays at a number of agar concentrations spanning the 
soft range (0.1-0.4%) should therefore be part of standard protocol when 
screening for chemotaxis. 



Acknowledgements 

We acknowledge work by Jessica Cameron in the embryonic stages of this 
research, assistance by Sarah Spragg, and discussions with Rosalind Allen, 
Julien Tailleur, Davide Marenduzzo and with Gary Darken. OAC and 
WCKP were funded by the EPSRC EP/E030173 and EP/D071070. GPF 
was funded by an MRC New Investigator grant G0501107. MEG was funded 
by the Royal Society. 



E. coli in soft agar 



14 



A Modelling run-and-tumble chemotaxis in agar 

Using a microscopic model of run-and-tumble dynamics in one dimension 
(see (|4d ) and references therein) it can be shown that the bacterial diffusivity, 
//, and the chemotactic parameter, x, are given by: 

2v'^ — a+ 
H = — ; v-y = V— (Al) 

a+ + a- ^ a+ + a- 

where are the mean tumble probabilities for bacteria moving up (+) 
and down (— ) the substrate gradient, and v is the average run speed. Note 
that for symmetric bias, a"*" = a~ = a, fj, = v'^/ad = fiQ in d dimensions. 
Extension to the asymmetric case for d > 1 is cumbersome and here we 
formally work only in d = 1. (By writing the final results in terms of the 
correct d-dependence is, however, recovered in the symmetric limit). 

We connect the above expressions to the chemotactic response by mod- 
ifying previous work (fisi ? ) to account for the effect of agar. A bacterial 
run is an inhomogeneous Poisson process with rate 



at{t) = ao 



1- / dt'K{t-t')f^{x{t')) 



ao[l-A(t)] (A2) 



where the subscript t indicates tumbles and, as in the main text, ao is the 
tumble rate in the absence of bias and f^. is a function related to substrate 
concentration at position x via = s{x)/{s{x) + k^). The function K{t) 
is the bilobed chemotactic response function which has been measured for 



E. coli (4^), and obeys Jq°° K{t)dt = 0. The linear expression above is valid 
for shallow substrate gradients, i.e. the bias |A(i)| <ti 1. Considering a run 
starting at t = 0, in the absence of agar the probability density for a tumble 

occurring in the interval [t, t + dt] is given by at{t) exp ^— dt'atit')^ . We 
argue that since bacterial collisions with the agar can also be considered a 
Poisson process, the same probability density describes the occurrence of 
tumbles in agar if the tumble rate at is replaced by an effective rate: 

ae(t;C) = at + aA (A3) 

which comprises (independent) contributions from at = at{t; C), the tumble 
rate due to the intrinsic bacterial dynamics (modulated by any chemotac- 
tic response) and a a = a^iC) an additional collision rate with the agar 
(which also randomises swimming direction). Then the mean run duration 
for bacteria in agar (or other porous media) is given by 

T{C) = (I dttae{t;C)exp(- [ dt' ae{t' ; C))\ (A4) 

XJo \ Jo ) I paths 
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where (. . ■) paths denotes an average over all possible bacterial swimming 
paths (the suffix will hereafter be assumed), since the nonlocal contribution 
at to Oe is path dependent. Then, changing variables in the memory integral 
(jMj) by defining u = t-t', substituting (jA3]) into (jAi]) and recalhng \ A{t) \ <C 
1, we have: 

T{C) « + ao [ dte-"(^)* / / dt' [ duK{u)f^{x{t' - u))) (A5) 

where we have defined the unbiased tumble rate in agar 

a(C7) = ao + aA(C). (A6) 

The concentration function is related to imposed gradients by a Taylor 
expansion: 

- It)) « x(i - n)V/x + const. (A7) 

Recalling that that K{u) integrates to zero, the constant term does not 
contribute to integral (jASp . Thus, following a trick introduced by De-Gennes 
(0), 

we consider a 'single delay' response function of the form K(u) = 
AS{u — 9), so ()A5P becomes 

T{C) « ^ + AVf^ao (ite-"(^)* l^j'^ dt'x{t' -6)^. (A8) 

Next, again following (fisl ) (ignoring persistence and rotational diffusion, see 
(|4ll)) we notice that for times t — 9<0 before the start of a run, the position 
x{t — 6) is on average not correlated to the bacterial velocity along the run. 
On the other hand for t — > 0, we can write x{t — 9) = ^v{t — 9), where 
±v is the run speed up or down a gradient. (|A8P then becomes 



1 f °° 1 

T^{C)^^^±v |V/^| aoA dte-^^'^'-it - 9f (A9) 



and, after integrating by parts 



. So finally, for a general a distribution of delay times K{6) we have 

-| f'OO 
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Now we identify = so that we can use (|ATT]) and dAl]) to find, to 

leading order in |V/^|, 



QO 



d9K{9)e 



(A12) 



a{Cy " a{Cf 

Or, since the chemotactic sensitivity parameter x is defined by = xiC)Vf^ 



1 + 



; x{C) = - 



1 + 



«a(C) 



deK{9)e 



-a{C)6 



(A13) 

where we have expanded the agar tumble rate defined in ()A6p . Equations 
()A13P are the bacterial transport parameters in agar accounting for a chemo- 
tactic response nonlocal in time. In the absence of agar [C — > 0), the ex- 
perimentally measured values of the bacterial transport parameters are 
and XO; the in liquido diffusivity and chemotactic parameter. In this limit 
()A13P become the expressions derived by de Gennes (fisl ) 



/xo; xiC^O) = — d0/^(0)e""°^ = xo. (A14) 
ao Jo 

Using ([Mill , we can rewrite (jAlSp as 



^{C ^ 0) = — 
ao 



/i(C) = Ho 
where 



1 + 



ao 



xiC) = xo 



1 + 



ao 



aA{C) 



ao 



ao 



deK{9)e- 



ao9 



(A15) 



(A16) 



To solve the model presented in the main text, we require an explicit expres- 
sion of K{t) to evaluate (jAlGh . and thus (jAlSh . We use a recently proposed 
fit to the experimentally measured impulse response of E. coli (j3c 
write: 

/ 1 2 2 

1 - ^0 aot + -aot 



and 



(A17) 



where ao is the base tumble rate, Ao is a dimensionless constant and A'^o > 
is a normalisation constant whose value is unimportant, as it cancels out in 
the expression for I^. 

To facilitate the discussion of our results, we also evaluate two limiting 
expressions for the concentration dependence of the chemotactic parameter 
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in (jAlSp . For very low concentrations bacterial collisions with the agar are 
rare, aA(C) ^ 1 (q^(C) ~ oo); so, expanding to first order, (|A16p becomes 
1 - KUA/ao, where k = f^^ deK{e)e-'^°'^aoe / deK{e)e-'^o'^ . For 
large agar concentrations, on the other hand, collisions with the agar are 
frequent and confuse the chemotactic response. The effective tumble rate 
is so large compared to the natural one, a{C) ^ ckq, that K{6) can be 
approximated by K{0) in the numerator of ()A16p . where the integrand falls 
rapidly to zero for 6 > l/a(C). In this case /-^(a^i/ao 3> 1) ~ /3[1 + 
aA(C)/ao]"^ where /3 = K{0)/{ao deK{6)e-''o'^). We can then write 
asymptotic expressions for the chemotactic parameter: 



X(C) 




,qa(C) 



if a{C) « ao 
if a{C) 3> ao . 



(A18) 



If, as in the main text, the values = 0.5 and qq = 1 are used in (|A17p . 
then K = 1/10 and (3 = 16/5. 
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Figure Legends 

Figure [H 

(a) Early (top row) and advanced (bottom) stages of the migration of E. 
coli AB1157 populations through LB agar of concentration C = 0.15-0.35%, 
as labelled. Shown are circular views (65 mm diameter) from minimally 
processed images (see methods) of 100 mm diameter petris filled with 10 
mm thick agar, (b-e) Azimuthally averaged radial intensity profiles from 
the images (see text). The time since inoculation in hours is indicated 
throughout. 



Figure [2], 

Bacterial populations for C = 0.4-0.5%. Colonies (65 mm views) are initially 
circular (top row), but quickly develop blebbing instabilities (bottom row). 
Images were minimally processed as for Fig. [1^ (see methods). 



Figure [3l 

Colony radius, r, against the time, t, since growth inception (see text) for 
C = 0.15-0.4%, as shown. The inset shows a linear fit to the raw radial data 
for C = 0.3% against time since inoculation, tj. Similar fits for all other 
concentrations define the migration speed (slope) and the latency time Ati 
(intersection with the time axis). Error bars are at most the size of a data 
point. 



Figure [H 

Experimental migration front speed as a function of concentration in the 
range 0.15-0.35% together with a best fit to the data using our model. Also 
shown for best fit parameters are the model 'efficient' and 'confused' limits, 
and the prediction from gas kinetic models (j23l ). The inset shows the same 
data but including points for C = 0.4-0.5%, labelled differently to indicate a 
different mode of migration at these concentrations. The model breakdown 
in this region is evident. 



Figure [5l 

Theoretical predictions for the band profiles for the full model (left) and gas 
kinetic models (j2^ ) (right) in the same range probed in experiments: C = 
0.15-0.35%, as indicated. In the full model, as concentration is increased 
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the dynamics changes from chemotactic (sharp bands) to growth/diffusion 
dominated (broad bands). This gradual transition is quahtatively the same 
as observed in experiment (see Fig. [1]), and is not predicted by the gas 
kinetic model. 
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